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Abstract 

The cross section of the process e"'"e~ Z77 — > qq77 is measured with 215 pb~^ 
of data collected with the L3 detector during the final LEP run at centre-of-mass 
energies around 205 GeV and 207 GeV. No deviation from the Standard Model 
expectation is observed. The full data sample of 713 pb~^, collected above the Z 
resonance, is used to constrain the coefficients of anomalous quartic gauge boson 
couplings to: 

-0.02 GeV~2 < ao/A^ < 0.03 GeV"^ and 
-0.07 GeV~2 < aJA^ < 0.05 GeV-^ 

at 95% confidence level. 



Submitted to Phys. Lett. B 



Introduction 



High energy e"'"e~ collisions offer a unique environment to unveil the structure of the couplings 
between gauge bosons. Extensive studies of boson pair-production are performed to probe triple 
vertices of neutral and charged bosons. Results were recently reported on the investigation of 
triple boson production through the reactions e+e~ W+W"7 [HE] and e+e" Z77 
These processes give access to possible anomalous Quartic Gauge boson Couplings (QGCs). 

Figures^— c display three of the six Standard Model diagrams that describe the e"'"e~ — 
Z77 process with the radiation of photons from the incoming electrons. This process is studied 
exploiting the high branching fraction of the Z boson decay into hadrons. The e"'"e~ — > Z77 — >■ 
qq77 signal is defined pO by phase-space requirements on the energies and angles 9^ of the 
two photons, on the propagator mass and on the angle between each photon and the 
nearest quark: 

E^>5GeV, I cos^^l < 0.97, jv^ - mz| < 2rz, and cos ^^.q < 0.98, (1) 

where mz and Fz are the Z boson mass and width. Events with hadrons and initial state 
photons falling outside the signal definition cuts are referred to as "non-resonant" background. 

A single initial state radiation photon can also lower the effective centre-of-mass energy of 
the e^e~ collision to around mz- This photon can be mistaken for the most energetic photon 
of the signal and two sources can then mimic the least energetic photon: the direct radiation of 
photons from the quarks, or photons originating from hadronic decays, misidentified electrons 
or unresolved vr^'s. These background processes are depicted in FiguresUJi andU^, respectively. 

In the Standard Model, the Z77 production via QGCs is forbidden at tree level. Possible 
contributions of anomalous QGCs, through the diagram sketched in Figure dF, are described 
by two terms of dimension-six in an effective Lagrangian 

TCfy — * — * 

where a is the fine structure constant, F^^, is the photon field and Wa^ is the weak boson 
field. The parameters oq and Oc describe the strength of the QGCs and A represents the scale 
of the New Physics responsible for these anomalous contributions. In the Standard Model, 
do = CLc = 0. Experimental limits on QGCs were derived from studies of the e"'"e~ —* W^W~7 
process However, the Oq and Oc couplings might be different in the e"'"e^ Z77 case. 

Alternative parametrisations can be found in References and ^ Indirect bounds on QGCs 
were extracted in Reference IHI using Z pole data. 

Data Analysis 

Reference m describes the analysis of the e"'"e~ — > Z77 qq77 process with 497.6 pb~^ of data 
collected by the L3 detector [12] at LEP at centre-of-mass energies, -y/i, between 130 and 202 
GeV. This Letter details the equivalent findings from the final LEP run, when the machine was 
operated at y/s = 200 — 209 GeV. These data are grouped in two energy bins around average 
y/s values of 204.8 GeV and 206.6 GeV, respectively corresponding to integrated luminosities 
of 77.4 pb"^ and 137.9 pb'^ 
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The signal and the "non-resonant" background are described with the KK2f Monte Carlo 
program ^Tj, which takes into account the interference of diagrams with initial and final state 
photons. It is interfaced with the JETSET ^2] program for the simulation of hadronisation. 

Other backgrounds are generated with the Monte Carlo programs PYTHIA |12| (e"'"e~ — *■ 
Ze+e- and e+e" ^ ZZ), KORALZ (e+e" ^ r+r-(7)), PHOJET [Hj (e+e" ^ e+e" 
hadrons) and KORALW [11] for W"'"W~ production except for the ez/eQq' final states, gener- 
ated with EXCALIBUR (TH]. The L3 detector response is simulated using the GEANT [TJj 
and GHEISHA jTHj programs, which model the effects of energy loss, multiple scattering and 
showering in the detector. Time dependent detector inefficiencies, as monitored during data 
taking periods, are also simulated 

Candidates for the e"'"e~ — » Z77 — ^ qq77 process are longitudinally and transversely bal- 
anced hadronic events with two isolated photons with reconstructed energy above 5 GeV, de- 
tected in a polar angle range | cos6'| < 0.97. The invariant mass of the reconstructed hadronic 
system, Mqq, is required to be consistent with mz- 74 GeV < Mqq < 111 GeV. 

The main background after these requirements is due to the "non-resonant" production of 
two photons and a hadronic system. The relativistic velocity jSz = pz/Ez of the Z candidate is 
calculated from the kinematics of the observed photons, assuming its mass to be mz- As shown 
in Figure 2a, f3z is larger for part of these background events than for the signal. Requiring 
Pz < 0.73 rejects half of this background. 

Events with a single initial state radiation photon, such as those shown in Figure Id and 
Figure le, are rejected by an upper bound on the energy E^i of the most energetic photon. This 
cut is chosen as E^i < 79.9 GeV at = 204.8 GeV and E^i < 80.6 GeV at ^/i = 205.6 GeV. 
A lower bound of 17° on the angle u between the direction of the least energetic photon and 
that of the closest jet is also imposed. Data and Monte Carlo distributions of these selection 
variables are presented in Figure |21 Good agreement is observed. 

Table 1 lists the signal efficiencies and the numbers of events selected in the data and Monte 
Carlo samples. A signal purity around 75% is obtained. The dominant background consists 
of hadronic events with photons. Half of these are "non-resonant" events, the other half being 
events with final state radiation or fake photons. 

Cross Section Measurement 

A clear Z signal is observed in the spectrum of the recoil mass to the two photons, as presented 
in Figure |HK- The e^e~ — > Z77 —>■ qq77 cross section, a, is determined in the kinematical 
region defined by Equation (1) at each average ^/s by a fit to the recoil mass spectrum. The 
background predictions and the signal shape are fixed, while the signal normalisation is fitted. 
The results ar 

a(204.8 GeV) = OMtoH ± 0.03 pb ((Tsm = 0.287 ± 0.003 pb) 
or(206.6 GeV) = 0.25l[J:[J^ ± 0.03 pb ((Tsm = 0.281 ± 0.003 pb). 

Here and below, the first quoted uncertainties are statistical and the second ones systematic. 
The systematic uncertainties on the cross section measurement are of the order of 10% 
mainly due to the limited Monte Carlo statistics and the uncertainty on the energy scale of the 
detector. 

The cross section is also measured in the more restrictive phase space defined by tightening the bounds on 6^ 
and d^q to I cos 6'-y I < 0.95 and cos6'-yq < 0.9. For the fuU 215 pb~^ at the combined average ^/s of 205.9 GeV, the 
resuh is: cr(205.9 GeV) = 0.18 ± 0.06 ± 0.02 pb, with a Standard Model expectation of ctsm = 0.172 ± 0.003 pb. 
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The measurements are in good agreement with the theoretical predictions, usm, as calculated 
with the KK2f Monte Carlo program. The uncertainty on the predictions (1.5%) is the quadratic 
sum of the theory uncertainty ^I] and the statistical uncertainty of the Monte Carlo sample 
used for the calculation. These results and those obtained at lower centre-of-mass energies jl] 
are compared in Figure El to the expected Standard Model cross section as a function of y/s. 

Figure Ob shows the recoil mass spectrum for the total data sample of 712.9 pb~^ collected 
at LEP above the Z resonance, comprising the data discussed in this Letter and those at lower 
centre-of-mass energies jl]. A fit to this spectrum determines the ratio Rzjj between all the 
observed data and the signal expectations as: 

R^^^ = ^ = 0.86 ±0.09 ±0.06, 

C"SM 

in agreement with the Standard Model. The correlation of systematic uncertainties between 
the different data samples amounts to 50% and is taken into account in the fit. 



^(GeV) 


e{%) 


Data 


Monte Carlo 


Ns 




j\qOtheT 


204.8 


51 


17 


14.7 ± 0.5 


11.3 ± 0.5 


3.09 ± 0.02 


0.31 ± 0.03 


206.6 


50 


23 


24.7 ± 0.5 


19.5 ± 0.5 


4.53 ± 0.04 


0.67 ± 0.03 



Table 1: Results of the e"'"e~ — >• Z77 — > qq77 selection. The signal efficiencies, e, are 
given, together with the observed and expected numbers of events. Expectations for 
signal, Ns, hadronic processes with photons, A'^^''^, and other backgrounds, Nether ^ 
are listed. Uncertainties are due to Monte Carlo statistics. 



Constraints on Quartic Gauge Boson Couplings 

Anomalous values of QGCs would manifest themselves as deviations in the total e"'"e~ — > Z77 
cross section as a function of y/s, as presented in Figure |3] A harder energy spectrum for the 
least energetic photon P constitutes a further powerful experimental signature, as shown in 
Figure El for the full data sample collected at ^/s = 130 — 209 GeV. QGC predictions for the 
cross section and this spectrum are obtained by reweighting the Standard Model signal Monte 
Carlo events. A modified version of the WRAP [19^ Monte Carlo program, that includes the 
QGC matrix element, is used. 

The energy spectra of the least energetic photon are fitted for the two y/s values discussed 
in this Letter and the eight values of ^/s of Reference HI Each of the two parameters describing 
the QGCs is left free in turn, the other being fixed to zero. The fits yield the 68% confidence 
level results: 

ao/A^ = 0.00i°:°2 GeV-2 and aJA^ = 0.031°;°^ GeV'^ , 

in agreement with the expected Standard Model values of zero. A simultaneous fit to both 
parameters yields the 95% confidence level limits: 

-0.02 GeV-2 < ao/A^ < 0.03 GeV^^ and - 0.07 GeV^^ < aJA"^ < 0.05 GeV'^ , 

as shown in Figure IHl A correlation coefficient of —16% is observed. Experimental systematic 
uncertainties as well as those on the Standard Model e"'"e~ —>■ Z77 —>■ qqyy cross section 
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predictions are taken into account in the fit. These resuhs supersede those previously obtained 
at lower ^/s 4J, as they are based on the full data sample and an improved modelling of QGC 
effects. 

In conclusion, the e~'"e~ — >■ Z77 — »• qq77 process is found to be well described by the 
Standard Model predictions JT], with no evidence for anomalous values of QGCs. 
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:+ 7 7 q e+ q 




Figure 1: Representative diagrams of, a)-c), the Standard Model contribution to 
the e+e~ — > Z77 signal and the "non- resonant" background, d), the background 
from direct radiation of a photon from the quarks, e), the background from pho- 
tons, misidentified electrons or unresolved tt^'s originating from hadrons and, f), the 
anomalous QGC diagram. 
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Figure 2: Distributions of, a), the relativistic velocity j3z of the Z boson reconstructed 
from the measured photons, b), the invariant mass Mqq of the hadronic system, c), 
the scaled energy E^i/ ^ of the most energetic photon and, d), the angle u between 
the least energetic photon and the nearest jet. Data, signal and background Monte 
Carlo samples are shown. Monte Carlo predictions are normalised to the integrated 
luminosity of the data. The arrows show the positions of the final selection cuts. In 
each plot, cuts on all other variables have been applied. 
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Figure 3: Mass recoiling from photon pairs in data, signal and background Monte 
Carlo for, a), the data sample analysed in this Letter and, b), the total sample 
collected above the Z resonance. Monte Carlo predictions are normalised to the 
integrated luminosity of the data. 




Vs (GeV) 

Figure 4: The cross section of the process e"'"e~ — > Z77 qq77 as a function of 
y/s. The signal is defined by the phase-space cuts of Equation (1). The width of the 
band corresponds to the statistical and theoretical uncertainties of the predictions 
of the KK2f Monte Carlo. Dashed and dotted lines represent anomalous QGC 
predictions for oq/A^ = 0.05 GeV~^ and ac/A^ = 0.10 GeV"^, respectively. The 
inset presents three combined samples: 231. 6pb^^ at y/s = 182.7 — 188.7 GeV, 
232.9pb~^ at y/s = 191.6 — 201.7 GeV and the data described in this Letter. 
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Figure 5: Energy spectrum of the least energetic photon in data, signal and back- 
ground Monte Carlo. The full integrated luminosity at ^/s = 130 — 209 GeV is 
considered. Monte Carlo predictions are normalised to the integrated luminosity of 
the data. Examples of anomalous QGC predictions are also given. 
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Figure 6: Two dimensional confidence level contours for the fitted QGC parameters 
Oo/A^ and Oc/A^. The fit result is shown together with the Standard Model (SM) 
predictions. 
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